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ABSTRACT 

Context. SN 1987A in the Large Magellanic Cloud is close enough for a study of the very late 
time evolution of a supernova and its transition to a supernova remnant. Nearly two decades after 
explosion we are witnessing the supernova shock wave engaging the inner circumstellar ring, 
which had been fluorescing since being ionised by the soft X-ray flash from shock breakout. 
Aims. We follow the interaction of the supernova shock with the ring material. The spatially re- 
solved information provides us with insight into the individual shock regions around the ring. 
Methods. Near-infrared integral field spectroscopy observations with SINFONI/VLT of the SN- 
ring interaction is presented. SINFONFs adaptive optics supported integral field spectrograph 
spatially resolves the ring and the data thus we obtain a better spatial understanding of the spec- 
trum in different regions of the object. 

Results. With a dynamical map of the interacting ring we determine parameters for its geometry. 
Since most of the IR emission lines originate behind the shock front we obtain an indication of 
the radial velocity of the shocked material after deconvolving the geometry. The ring geometry is 
consistent with a circle and we also derive a new, independent measurement of the systemic ring, 
and presumably also supernova, velocity. We find from the spatial distributions of the flux in the 
different emission lines the degree of cooling in the shocked material and follows the increases 
observed in the radio and X-rays. Emission from the ejecta is detected only in the strongest [Fe n] 
lines. 
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1. Introduction 

With the 8m telescopes in the Southern Hemisphere SN 1987 A will remain observable through- 
out its complete evolution from supernova to supernova remnant. While t he emission from the 



supernova itself has been steadily fading (e.g. lLeibundgut & Suntzefl 



2003) the inner circumstel- 



lar ring has remained bright. The collision of the supernova shock with the circumstellar ring has 

j II — 1 1 r~i 

been anticipated (Luo et al. 1994; Borkowski et al. 1997) and was first indicated by radio observa- 



tions (Ball etal. 2001 



Manchester et al. 2002). At about the same time the first hot spots appeared 



just inside the circumstellar ring in 1996 showing that the supernova shock had reached what are 



Based on observations collected at the European Southern Observatory, Chile. 
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considered to be inwar d protrusions of the ring dSonneborn et al 



1998 



Michael et al 



2002; 



Pun et al 



Haberl et al 



20021) . The X-ray emission started to increase as well dPark et al 



1998 



2004 



2000 



2006 



20061) . The detection of coronal lines in the opti cal spectra from 2002 in dicates that 



the ring material is sho cked to a temperature of :> 2 x 10 6 K dGroningsson et al, 



with the X-ray spectra dZhekov et al 



2005 



2006). 



2006), consistent 



The physics of the interaction is highly complex because of the nature of the circumstellar 



medium of SN 1987 A dChevalier & Dwarka das 1995 



Lundqvist 1999). 



Here we present spatially resolved spectroscopy of the inner circumstellar ring from science 
verification observations with SINFONI on the VLT. Integral field spectroscopy yields both imag- 
ing and spectroscopy, which allowed us to trace the shock interaction in spatially separated emis- 
sion sites around the ring. The approaching and receding sides can be separated very clearly in 
several emission lines and we are able to study the ring dynamics. 

The optical and IR provide, together with the UV and X-rays, the most important diagnostics 
of the temperature, density, and ionisation of the shocked gas. Therefore, the data presented here 
should be seen as a complemen t to the high resoluti on optical and UV spectr oscopy presented in 



Groningssonetal 
2006). 



(2006, 



2007) and 



Pun et al 



(2002), as well as the X-rays by 



Zhekov et al 



(2005 



2. Observations and Calibration 

We use SINFONI Science Verification observations of SN 1987A carried out on 28 and 29 



November 2004 (day 6488, preliminary results were published in Gilles sen et al.H 2005). Integral- 
field spectroscopy was obtained in J (1.1 1-1 .35 yum), H ( 1 .49- 1 .79 fim) and K ( 1 .99-2.38 fim) bands. 
The observations were AO supported using Star 3 as the referenc e source. We h ave not deconvolved 



the images, since deconvolving extended sources is unrealiable (Davies 2007). The pre-slit optics 
for SINFONI allows to chose different widths of the image slicing, which then leads to different 
fields of view on the sky. SN 1987A was observed with two different spatial resolution settings. 
The primary data in J, H, and K were obtained with a spatial resolution of 250 mas/spaxel and thus 
with a field of view of 8"x8". In J we also obtained an observation with a spatial resolution of 100 
mas/spaxel and a field of view of 3"x3" (denoted as J100 in the following). All single integrations 
were 600 seconds and Table [T]displays the total integrations per setting. When several integrations 
were obtained the individual exposures were offset by sub-spaxel spacings to improve the spatial 
resolution. The exposures were separated by sky exposures of equal integration time in an ABBA 
sequence. 

The sizes of the spaxels for the observations are 250 mas x 125 mas (and 100 mas x 50 mas 
for J100). In the K-band images the final image quality is superior to the J and H band due to better 
seeing and sampling. This can been seen in FigureQ] which shows the images (collapsed cubes) for 
the different bands. In the final combined data cube the spaxels are square and have the sizes 125 
mas x 125 mas (50 mas x 50 mas for the J100). The ring is fully resolved in all bands. Using the 
stars in the FOV of the collapsed images we measure an image quality (values listed in Table [1). 
The image quality measurement combines seeing, the AO correction and any smearing/dithering 
in combining the final cube. The image quality is in particular important in estimating reasonable 
sizes of spectral extractions (apertures). 
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Fig. 1. The 8"x8" FOV in the J, H and K-band images. Lower right is the higher resolution J-band 
3"x3" FOV. 



The data have been reduced and recom bined into data cubes using the SINFONI pipeline 



( Schreiber et al 



2004 



Modigliani et al 



2007). We crosschecked the wavelength solution with the 



atmospheric OH lines dRousselot et al.l 2000) and the resulting wavelength calibration is accurate 
to 8x10 ~ 5 /im, 7xl0~ 5 //m and 7.5x10 -5 /mi, in J, H and K, respectively. This corresponds to sys- 
tematic velocity uncertainties of 19 km s , 12 km s _1 , and 1 1 km s _1 , as li sted in Tabled 
We dereddened the spectra using t he galactic extinction law fr om 



Cardelli et al 



3.1, and Eg_y=0.16 dFitzpatrick & Walborn 



(1989), 



1990) for the color excess to- 



assuming R v 

wards SN 1987A, based o n E B - V =0.10 from the LMC and £ B _y=0.06 from the Milky Way 



dStavelev-Smith et al 



20031) . The differences between the LMC and galactic extinction law are 



negligible at low color excess in the near-IR. 

The standard stars were observed directly before or after the supernova and have been used to 
remove the telluric features. We used the standard stars to construct a sensitivity cur ve and derive 



the flu x calibration, based on their measured magnitudes from the 2MASS catalogue dCohen et al 



2003). In this way the whole data cube had the same flux calibration, and we could derive 
calibrated spectra for Stars 2 and 3 as well. T o test the absolute flu x calibration we compared the 
spectrum of Star 2 (a B0V star NW of the SN, Walborn et"alll993 ) with the blackbody curve of a 
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Table 1. Summary of observations. The rows designate the observation date (UT), the time since the 
explosion of SN 1987 A on 23 February 1987 A, the spaxel size of the observation, the wavelength 
range of the spectrograph setting, the spectral resolution, the total exposure time and the image 
quality as measured from the combined data cube on Star 2 and Star 3. We further list the standard 
star used for the flux calibration, and the uncertainty in the flux as determined from a comparison 
of Star 2 from the corresponding black body curve, the uncertainty in the wavelength calibration 
as measured from night sky lines and the corresponding uncertainty in velocity. 





J 


H 


K 


J100 


Date 


2004/11/28 


2004/1 1/29 


2004/11/29 


2004/11/28 


Day 


6488 


6489 


6489 


6488 


Spatial Res. 


250mas 


250mas 


250mas 


lOOmas 


Wavelength 


1.11-1.35 iim 


1.487-1.792 yum 


1.985-2.375 fim 


1.11-1.35 fim 


A/AA 


2000 


2000 


4500 


2000 


Exptime 


600s 


1200s 


1800s 


1200s 


Image quality 


0.62" 


0.46" 


0.48" 


0.32" 


Std. Star 


HD36531 


HD29303 


HD29303 


HD36531 


EFlux 


1 % 


13 % 


4 % 


4% 


S(A) 


8 xlO^/mi 


7 xlO^/mi 


7.5 xl0~ 5 fim 


8xlO~ 5 /im 


<S(v) 


19kms-' 


12kms-' 


11 kms- 1 


19 kms 1 



star with R=8.3 R Q and T e //=22500 K at the distance of 50 kpc. The derived flux differences are 
of the order of 10% and are listed as €pi ux in Table [TJ For the J 100 observations we calibrated the 
flux using Star 3, since Star 2 is not in this smaller FOV. As Star 3 is variable we used the observed 
values of Star 3 in the J250 to cross correlate the value. The total flux from Star 3 observed in J 100 
was 4% different from the total flux from Star 3 obtained from J250 using the same aperture in 
arcseconds. We combined the two J observations into one spectrum weigthing the spectra with the 
exposure time. The flux error estimate for this combined spectrum is of the order of 2%. We do 
not think that any emission contamination from the ring will influence the result beyond the 10% 
we find as the accuracy of the flux. 

3. Results 

Because the data format is 3-dimensional (spatial, spatial, and wavelength) it is possible to integrate 
spectra over different spatial regions. Generally, it is only correct to compare line fluxes with respect 
to other line fluxes, where the number of co-added spaxels, i.e. the aperture, is the same. In order to 
compare with previous data we summed up the circumstellar ring and the ejecta (which corresponds 
to 28 1 spaxels co-added for the 250 mas data cubes). We extracted the spectrum in a similar way for 
the J100 observations and scaled this spectrum to the J spectrum before we combined the spaxels 
into a single spectrum. 

We chose to scale the J 100, because the aperture for the J band is easily compared with the ones 
of the H and K bands, which were observed with the same spatial scale. This spectrum should be 
comparable with earlier spectra of SN 1987 A after a slit correction. We will refer to the spectrum 
extracted this way as the 'integrated spectrum' and it is displayed in Figure [2] We also show a 
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Fig. 2. The integrated spectrum of the ring and the ejecta. The sky spectrum is shown offset below 
the SN ring spectrum. Note that the flux of the sky spectrum has been scaled down to make the 
spectrum fit in the same figure as the object. The line identifications are given in Table|2]|4] 



sky spectrum offset below the 'integrated spectrum'. This sky spectrum should be seen as a help 
to estimate where the more noisy regions are located in wavelength. The intensity of the various 
skylines cannot be directly compared with the integrated spectrum due to differences in aperture. 

A very weak continuum flux can be discerned in the spectrum at (2 - 3) x 10 _17 ergs _1 cni _2 A _1 . 
The consistency of this continuum in all three bands makes us believe that it is real. The origin of 
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Table 2. Emission lines of the integrated spectrum (ejecta and circumstellar ring combined corre- 
sponding to the central 281 spaxels) in the J band. The spectral resolution is 150 km s in the J 
band. Errors for the values listed in the table are given either in the text or in TableQ] 



Aobs 




O" obs 


(Telecom* 


Flux 


Identification 


■lair 


(//m) 


(km s" 1 ) 


(kms- 1 ) 


(km s-') 


(10- 16 ergs-' cm- 2 ) 






1.253885 


408 


431 


404 


14 


*[Fen] a 6 D 3/2 -a 4 D 1/2 


1.2521 


1.25385 


261 


431 


404 


14 


*He I ls3s-ls4p(3S-3P) 


1.2528 


Blend? The line transition probability for [Fe n] does not support this large flux. Furthermore large cr deconv . 


1.25787 


277 


200 


132 


69 


f[Fen] a 6 D 9/2 -a 4 D 7/2 


1.2567 


1.27169 


315 


236 


182 


7 


f[Fen] a 6 D 1/2 -a 4 D 1/2 


1.2703 


1.27983 


245 


305 


266 


43 


t[Fen] a 6 D 3 / 2 -a 4 D 3/2 


1.2788 


1.27983 


310 


305 


266 


43 


*HeIls3d-ls5f(3D-3F) 


1.2785 


Blend? The line transition probability for [Fe n] does not support this large flux. Furthermore large cr deconv . 


1.28300 


277 


281 


238 


367 


fPa/3 


1.2818 


1.29562 


311 


190 


117 


18 


t[Fen]a 6 D 5/2 -a 4 D 5/2 


1.2943 


1.30651 


302 


253 


204 


12 


*Ni 


1.3052 


1.32178 


277 


166 


71 


23 


t[Fen] a 6 D 7/2 -a 4 D 7/2 


1.3206 


1.32901 


226 


339 


304 


10 


t[Fen] a 6 D 3/2 -a 4 D 5/2 


1.3278 



The lines marked with t have been observed previously for this object in Meikle et al. (1993). 
The lines marked with * are new for this object or have a new identification. 



this continuum is not clear, but it could come from large scale emission in the region. If it does 
indeed emerge from the ring system itself then it could be a weak H/He continuum. Most likely 
the continuum is a combination of several emission origins, and it is at the moment not possible to 
distinguish between them. 

The detected lines from the spectrum are listed in Tables[2][3]and[4] The weakest lines we have 
been able to securely identify are those with a flux above 3xl0~ 16 erg s _1 crrT 2 . In the tables A ODS 
denotes the observed wavelength, X air is rest wavelength in air for the line transition listed under 
'Identification'. The heliocentric velocity of a given line is listed in column v i, s and calculated as 
the velocity shift between the observed wavelength and the rest wavelength for the identification 
and then corrected for heliocentric velocity. cr ois is the observed FWHM of the line and o~d econv 
refers to the FWHM for the line corrected for the instrument resolution {crdeconv - -^cr^ 5 - cr 2 instr ). 
The integrated flux for the line is measured by fitting a Gaussian to the line shapes and listed 
in the column labelled 'Flux' in the tables. We chose to ignore unresolved narr ow fines, which 
we co uld identify as residue sky subtraction lines using the OH line catalogue of 



Rousselot et al 



(2000). Since we have a 3D data format we could investigate the spatial extent of each line in 
the integrated spectrum and check whether the line originated at a particular point in the object or 
whether it has a coherent appearance across the ring/ejecta. This is a very effective way to check 
for artefacts, such as hot/bad pixels. All observed values are integrations over the whole ring. 

We have used the NIST Atomic Spectra Database Version 3.1.0. (Martin et al .120031) to identify 
the lines, checking the following elements and their ions: H, He, C, N, O, Na, Mg, Si, S, Ca, Fe, 
Co, and Ni. The ions included in the NIST database are up to ionisation XXIV for the case of iron. 
The other elements are complete to a similar level. The few lines which we have not yet identified 
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Table 3. Emission lines in the integrated spectrum (ejecta and circumstellar ring combined corre- 
sponding to 281 spaxels) in H. The spectral resolution is 150 km s in the H band. Errors for the 
values listed in the table are given either in the text or in Table [TJ 





Vobs 


O~ bs 




Flux 


Identification 




(jim) 


(kms- 1 ) 


(km s _1 ) 


(km s- 1 ) 


(10-' 6 erg s" 1 cm" 2 ) 




(jim) 


1.49825 


301 


320 


283 


2 


*HlBr25 


1.4967 


1.50986 


295 


278 


234 


5 


tHelSs'S^p'P 


1.5084 


1.51461 


253 


257 


209 


3 


*HlBr21 


1.5133 


1.52037 


232 


276 


232 


3 


*HlBr20 


1.5192 


1.52705 


197 


334 


298 


4 


*HlBrl9 


1.5261 


1.53494 


285 


254 


205 


19 


t[Fen]a 4 F 9/ 5-a 4 D 5 /2 


1.5335 


1.54512 


237 


272 


227 


5 


*HlBrl7 


1.5439 


1.55707 


273 


270 


224 


5 


*HlBrl6 


1.5556 


1.57142 


257 


305 


266 


7 


fHlBrl5 


1.5701 


1.58947 


273 


283 


240 


10 


tSii4s 1 P°-4p 1 P 


1.5880 




267 








fHlBrl4 


1.5881 


The transition probability for Brl4 alone fits with the flux. Sii identification tentative. 




1.60098 


281 


225 


168 


11 


t[Fen] a 4 F 7/ 2-a 4 D 3/ 2 


1.5995 


1.61232 


257 


279 


235 


12 


fHlBrl3 


1.6109 


1.64236 


298 


219 


160 


12 


fHlBrl2 


1.6407 


1.64438 


150 


3393 


3390 


47 


t[Fen]a 4 F 9/2 -a 4 D 7/2 


1.6435 


Broad component. See text for comments. 








1.64518 


296 


171 


82 


51 


t[Fen] a 4 F 9/2 -a 4 D 7/2 


1.6435 


This component sits 


on top of the broad component (previous line). 






1.66534 


282 


234 


180 


6 


t[Fen]a 4 F 5/2 -a 4 D 1/2 


1.6638 


1.67851 


290 


233 


178 


12 


t[Fen] a 4 F 7/2 -a 4 D 5/2 


1.6769 


1.68227 


286 


268 


222 


15 


fHlBrll 


1.6807 


1.70171 


259 


247 


196 


15 


tHei3p 3 P°-4d 3 D 


1.7002 


1.71289 


307 


263 


216 


3 


t[Feii]a 4 F 5/2 -a 4 D 3/2 


1.7111 


1.73769 


254 


311 


272 


21 


fHlBrlO 


1.7362 


1.74661 


286 


327 


291 


6 


t[Feii]a 4 F 3/2 -a 4 D 1/2 


1.7449 



The lines marked with t have been observed previously for this object. 



The lines marked with * are new for this object or have a new identification. 



did not correspond to any of the mentioned elements using the NIST database version 3.1.0, which 
is not necessarily complete for these elements and all their ions. The accuracies of the wavelengths 
in NIST are implied by the number of significant figures. 

For the [Fe n] and H i lines we h ave verified the identif i cation by checking line ratios using 
the line transition probabilities from iNussbaumer & Storevl (11988) and the NIST database. The 



broad component in the H-band (Table [3} is the [Fen] 1.644 jjm line originating in the ejecta. The 
flux from the narrow component sitting on top of the broad component alone agrees well with 
the fluxes observed from the 1.257 /urn and 1.321 /urn [Fen] lines, which come from the same 
upper level. However, if the broad component is the [Fen] 1.644 jum line then we should expect 
to see broad components from the lines [Fen] 1.257 fim and 1.321 /im as well. Given the S/N in 
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Table 4. Emission lines in the complete spectrum (ejecta and circumstellar ring combined corre- 
sponding to 281 spaxels) in K. The spectral resolution is 67 km s in the K band. Errors for the 
values listed in the table are given either in the text or in Table [TJ 





v o6! 






Flux 


Identification 


A a j r 


(/an) 


(km s-') 


(km s" 1 ) 


(kins" 1 ) 


(10- 16 ergs-' 


cm -2 ) 


(jim) 


2.00855 


275 


239 


229 


5 


t[Fen] a 4 P 1/2 -a 2 P 1/2 


2.0067 


2.01696 


272 


327 


317 


8 


* [Fen] a 2 G 9/2 -a 2 H 1/2 


2.0151 


2.04779 


259 


293 


285 


8 


t[Fen] a 4 P 5/ 2-a 2 P 3/ 2 


2.0460 


2.06003 


275 


291 


283 


108 


tHei2s'S-2p'P° 


2.0581 


2.11412 


279 


369 


363 


8 


*Hei Is3p-ls4s 


2.1121 


2.13462 


258 


267 


258 


3 


t[Fen] a 4 P 3/2 -a 2 P 3/2 


2.1328 


2.16331 




416 


411 


4 


* 




Which could be part of 


a broad component of Bry (not seen on the red side of Bry). 




2.16754 


277 


291 


283 


58 


fH I Bry 


2.1655 


2.19052 




575 


571 


3 


* 




2.20659 


255 


381 


375 


2 


*Nai2p 6 5p - 2p 6 9s 


2.2047 


2.20659 


-58 


381 


375 


2 


orfNar4i 3 5 -Ap 2 P° 


2.2070 


2.22026 




365 


359 


4 


* 




2.22593 


290 


243 


234 


5 


* [Fen] a 2 G 9 /2-a 2 H 11/2 


2.2238 


2.24580 


286 


534 


530 


4 


t[Fen] a 4 P 1/2 -a 2 P 3/2 


2.2436 


Blend? The line transition probability for [Fen] does not support this large flux. Furthermore large crvfeanw 




2.25525 




452 


447 


3 


* 




2.31030 


304 


403 


397 


3 


*Nin 3p 6 3d\ 3 F)4s - 3p 6 3d s ( 3 F)4s 


2.3079 


2.32381 


291 


284 


276 


5 


*Oi2s 2 2p 3 ( 2 D)3s - 2s 2 2p 3 ( 4 S)5f 


2.3215 



The lines marked with t have been observed previously for this object. 

The lines marked with * are new for this object or have a new identification. 

The lines 2.24-2.32 yum were observed in Meikle et al. (1993) as a blended unidentified feature. 



our J observations they are not discernible. The two lines are located in a noisy region of the J- 
band especially when compared with the H-band region around the 1.644 fim line. This point can 
hopefully be resolved with higher signal to noise data. 

When the identification transition in Tables |2]|4] i s preceded w i th a da gg er (f), this tr a nsitio n 



Meikle et al 



( 1993|) or 



Fassia et al 



(2002). 



has been observed previously and was identified by 
Generally, the ring lines are now brighter than for those observations, as expected if the shocked 
component dominates the narrow component from the unshocked ring. With the higher spectral 
resolution of SINFONI we can now resolve some previously blended lines. However, there remain 
blends for lines with a true separation of less than 150 km s _1 . We mark newly identified lines with 

an asterisk (*) in the tables. 

From HS T spectroscopic observations (Pun et alJr2002 ) and our ground-based high-resolution 



spectroscopy dGroningsson et al. 



2006, 



2007), we know that the ring-ejecta system at the moment 



contains three different velocity components. There is the narrow component (NC), which displays 
a nearly Gaussian velocity distribution, with cr ~ 10 km s _1 , arising from the unshocked circum- 
stellar ring. This is caused by fluorescence f rom the material recombining after being ionised by 



the supernova's UV flash at shock breakout dFransson & Lundqvist 



1989; 



Lundqvist & Fransson 
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Fig. 3. North and South extractions from SINFONI J 1 00 of the [Fe n] 1 .257 fim, where the SINFONI 
data have been spatially smoothed to the image quality of the UVES observations, and extracted 
along a position angle (PA=30°) of the UVES observations. The red/dashed curve are North/South 
extractions from UVES of the [Fen]7155 A line reduced to the spectral resolution of SINFONI 
scaled up a factor 100. 



19911) . The intermediate component (IC), with cr ~ 200 - 300 km s originates from the shocked 
ring material behind the passing shock fron t. The brightening has been observed in HST imaging 



(Michaeletaiil998 



2000; 



Pun et al.l l2002) and arises from this material, where the ejecta collides 



with the inner protrusions of the ring. In our observations the NC and IC cannot be distinguished 
and appear as a single marginally resolved component. The SN ejecta component wit h an expan 



sion velocity of ~ 1500 km s 1 comes from the interaction with the reverse shock dSmifh et al 



2005 



Heng et al. 2006) not evident in our data. We also observe a broad component (BC) which 



will be discussed further in section I3T1 

Because we do not have the spectral resolution to separate the narrow from the intermediate- 
velocity component (cf. Tab. [TJ we expect some of the intermediate lines we observed to be a 
convolution of the narrow component with the intermediate component. In addition we have inte- 
grated spatially over the whole ring and thus the lines are also broadened by the velocity disper- 
sion of the ring. Figure [3] shows a comparison of [Fen] emission lines observed with UVES and 
SINFONI, where the UVES data has been scaled up by a factor of 100. The SINFONI observations 
have been convolved to the seeing of the UVES observations (0.5")- We extracted two spectra, one 
North and one South on the circumstellar ring emulating the 0.8" slit (PA=30°), which was the 
setup of the SN 1987 A observations with UVES. The UVES spectra from the North and South part 
were convolved with the spectral resolution of SINFONI (A/AA -2000 in the J-band) and shown 
in the figure as the dashed (red) curve. Taking into account that the wavelength accuracy in the 
J-band is + 19 km/s then we have a good agreement for the peak velocity between the two observa- 
tions. A detailed compar ison of the individual elements with the high resolution optical spectrum 
( Groningsson et al]|2007 ) will further help in the separation of the different lines. 




Fig. 4. The [Fen] A 1.644 /j.m line from a single spaxel centered on the ejecta at the center of the 
circumstellar ring (aperture=l spaxel). The intermediate/narrow line on top is scatter originating in 
the ring. The red/dashed curve is the extraction at 25° (aperture area A-n * (0.375") 2 ) scaled down 
with a factor of 10. 



3.1. The Broad Component 
3.1.1. Central Emission 

The only line, for which we clearly can separate the broad ejecta component from the intermediate- 
velocity component in the integrated spectrum is the [Fen] Al. 644 fim line. Figure |4] shows this 
line from a single spaxel centered on the ejecta in the center of the circumstellar ring(full curve), 
for comparison we have plotted an extraction from Spot 1 (dashed line). The emission from Spot 1 
has the intermediate BiT2 line which is not present in the ejecta emission, furthermore there seems 
to be more blue emission in Spot 1 than in the ejecta. 

For the integrated spectrum when we fit a Gaussian to both the broad part and the intermediate 
part of this line we find a small offset in the central wavelength, so that the broad line seems 
less red-shifte d than the intermed iate line, with respect to the rest frame of the system {v$n = 
286.5 kms" 1 , Meaburn et alll995 ). The two components are individually listed in Table[3] For the 
broad component we derived a velocity shift of 150 kms~' and a ov cram , = 3400 km s _1 . The bluer 
central wavelength for the ejecta comp onent cou l d ind icate that the red part of the line is obscured 
by dust, as observed in the optical by 



Lucy et al 



(1199 ll) 



This is also an indication on the faintness of the ejecta 17.8 years after the explosion. The ejecta 



spect 
2004, 



rum is now dominated by energy input from Ti dFransson & Kozma 
2008) . 



2002; 



Bouchet et al 



3.1 .2. Ring Emission 

We assigned an azimuthal angle for each position on the ring, defining North as zero degrees, and 
increasing the angle through East. Figure [5] shows the size of the extractions and their correspond- 
ing positions and azimuthal angles. Figure [6] shows the extracted spectrum of [Fen]1.644//m for 
different positions on the circumstellar ring and the zoom of three of the extractions. It is clear from 
this figure that the velocity shift of the [Fen] 1.644 /urn line's narrow component is not changing 
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radius 0.375" 




Fig. 5. The size of the extractions and the corresponding azimuthal positions in the H band. 



much for different azimuthal angles (different positions on the ring). Around Spot 1 and on the 
eastern side in general a broad component of the line emerges (~ 8000 km/s FWZI). The shape of 
the broad component (BC) changes even on the eastern side, with a slightly more luminous blue 
wing at 25 degrees, and a slightly more luminous red side at 99 degrees. However the ratio of blue 
to red flux in the broad component changes with azimuthal direction away from Spot 1 (see Fig. 
01. Therefore we conclude that the emission is not only scattered emission from Spot 1. Given 
the accuracy of this data se t we cannot defin itively exclude contamination. Following the thorough 
discussion of the shocks in 



Pun et al 



(2002) we think that the broad emission in the ring emerges 
from the reverse shock, which would be bluer in the North part of the ring due to the inclination of 
the ring with respect to the line of sight, which is consistent with our observations. Again, a deeper, 
high signal to noise observation would solve this question. 

The broad line is only pronounced on the Eastern side of the ring. This East - West asymmetry 



Park et al. 



(20061) and 



Manchester et al 



has al so been observed in X-rays and radio most recently in 
d2002l) . respectively. The X-ray and radio emission is thought to emerge from the same region as 
the IR behind the shock front as the material cools down. We see also an East- West asymmetry 
with the intermediate lines, in that the Eastern side is brighter than the Western (see Fig. [T0b. 

Whether the East- West asymmetry is due to asymmetric ejecta outflow or differences in the 
density of the circumstellar matter (protrusions) is still an open question, which integral field spec- 
troscopy observations of the reverse shock would shed light on. Since the hot spots primarily first 
appeared on the Eastern side we could also be witnessing spot evolution. This is a question which 
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Fig. 6. Extraction of spectra from the reduced data cube in the H band. The spectra are plotted with 
an offset. The faint emission on the blue side of [Fen]1.644/mi is Brl2. 



will be addressed in an upcoming paper which presents a higher spatial resolution data set for the 
epoch 2005. 

3.2. The Circumstellar Ring 

The emission lines from the circumstellar ring are a complex combination of different emission 
sites. With the shock moving into the ring material an outward acceleration of the shocked material 
is expected, and hence a difference in the p rojec ted radial veloci ty shift around the ring should 
be observed. However, 



Michael et al. 



(2000) and 



Pain et al 



(2002) showed that there are oblique 



shocks along the protrutions and the velocity field is a combination of material accelerated not 
only radially forward. The SINFONI data cube can be used to explore the velocity structure of the 
ring. We assigned an azimuthal angle for each position on the ring, defining North as zero degrees, 
and increasing the angle through East. The resolution in angle is then set by the spatial resolution. 
For each angle we coadded a spectrum with the angular position as the center and the radius of 
the order of the image quality. In this sampling two neighbouring extracted spectra therefore have 
many spaxels in common and are correlated (cf. Fig. |7). The radius of the extraction circle was 
3 spaxels for the 250 mas resolution (radius of 0.375") and 5 spaxels for the 100 mas resolution 
(radius of 0.25"). We have here chosen the radius to be of the order of the image quality listed 
in Table [T] for the different observations. The left panel of Figure [7] shows the extracted spectrum 
(J 100-band between 1.25^m-1.29^m ) for different positions on the circumstellar ring and the right 
panel shows the sizes of the extractions and their corresponding positions and azimuthal angles. It 
is clear that the peak of [Fe n] 1.257 jjm and the peak of Pa/? change for the different extractions. 

Figure[8]shows the radial velocity shift for the Pa/3, Bry, Hei, and [Fen]1.257/mi lines, plotted 
vs. the azimuthal angle. We have here chosen the extraction radii for the J100 extractions to be 
similar to that of K250 (radius=7 spaxels ~ 0.35"). The error bars in velocity are calculated from 
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271" 
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Fig. 7. Left: Extraction of spectra from the reduced data cube close to the [Fe n] 1 .257 yum and Pa/? 
lines in the J 100 band. The spectra are plotted with an offset. Right: The size of the extractions and 
the corresponding azimutal positions in the J 100 band. The vertical lines through the spectra mark 
the positions of the systemic velocity of 286.5 km/s for the different lines. 

the uncertainty in the determination of the wavelength and the error bars in angle display the size 
of the extraction radius for the spectra. The error bars in angle are plotted conservatively to indicate 
the extent of contamination between adjacent spaxels. 

Only extractions with a distance of more than twice the integration radius are truly independent 
of each other. For J, H, and K the radius was 3 spaxels, which corresponds to a separation of 6 
spaxels (0.75"), and 14 spaxels for J100 (0.7"). 

The emission lines in Fig. [8] are all dominated by the intermediate component, coming from 
the shocked ring material. This is clear from the fact that the velocity of the lines differs by a large 
factor from the expansion velocity of the narrow lines, ~ 10 km s . It is evident from this figure 
that the line centers of the northern part of the ring are blue shifted, while the southern part is red 
shifted. This is as e xpected for an e xpan ding circular ring with the Northern part tilted towards us. 



From HST imaging iPanagia et al 



( 1991) find an inclination angle i=42.°8 + 2.°6 with respect to 



the plane of the sky (see also lSugerman et all2005l) . The theoretical line of sight velocity, v, , for a 
circular ring, uniformly expanding with velocity v exp , and tilted with respect to the plane of the sky 
with the angle i, is 



v r (i//) = v eX p sm i cos OA + <P) + v 



SN 



(1) 



Here ifr is the azimuthal angle of the ring segment, defined with if/ — for North. The phase shift, 
<p, is introduced to account for the offset of the major axis of the ring from the East- West direction 
and can be interpreted as a rotation of the ring out of the plane of the sky. 

The observed velocity shift at a given angle is a complex convolution of the individual shock 
velocities at this point and th e geometry. A simple model of the line profile of an individual spot 



was given in 



Pun et al 



( 2002). The key point of this model is that the shock velocity is likely to vary 
along the surface of the hot spot depending on the angle between the surface and the impacting blast 
wave. The highest velocities are therefore expected for a head on collision, while the tangential 
impacts result in lower shock velocities. Furthermore, the gas velocity, in the reference system of 
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Fig. 8. Radial velocity shift for different angular positions around the ring for the strongest emission 
lines. The angle increases from North to East. The location of Spot 1 and Star 4 are indicated. The 
J 100 data has been sampled with 7 px radius in order to have a similar extraction radius as the 250 
obs. The errorbars are indicated in the top right corner. 



the shock, decreases behind the shock as the gas cools. For the observer at rest, the hot gas velocity 
behind the shock is 3VJ4, while that of the radiatively cooled gas is V s . The peak velocity at a 
given angle therefore depends on the shock velocity which dominates the contribution to a given 
line. It is therefore important to realise that the peak velocity we measure is only a weighted average 
of the emissivity of a given line. This is then modified by the angle of the gas velocity to the line 
of sight. 

This can clearly be seen in the VLT/UVES observations, which provide lin e profiles with much 



higher spectral, but lower spatial resolution dGroningsson et al.l 1200 6, 2007). The line profile of 



e.g., the Hydrogen Balmer and He I lines from the northern part of the ring, including Spot 1, has 
a peak at ~ 200 km s -1 in November 2005. This agrees well with that measured at the position of 
Spot 1 in our observations, ~ 223 + 7 km s _1 (combining Pa/?, Bry and Hel 2.058 yum). The Ha 
line in the UVES observations, however, extends out to ~ -450 km s , seen only as a faint wing. 
While the UVES observations provide us with a good representation of the emission as function 
of velocity, the SINFONI observations measures the variation of the average velocity for each line 
along the ring. This illustrates well the complementarity of the SINFONI and UVES observations. 

It is clear from Fig.|8]that the measured velocities do not follow the simple cosine as predicted 
by the Eq. Q. 

A simplistic way to derive the bulk velocity of the material, i.e. the center of rest for the ring, 
would be to determine the 'nodal' points {ijj = 90° and \j/ = 270°), where the expansion velocity is 
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orthogonal to the line of sight. Inspection of Fig.[8]immediately shows that the two velocities are not 
identical at these angles for the various lines. We find that for Pa/3, Bry and He i v r (90°) = (300 + 7) 
km s _1 and v r (270°) = (275 + 3) km s . For a perfect circular ring tilted out of the plane of the 
sky the radial velocity should vanish at points exactly 180° apart and these points indicate the bulk 
velocity of the material. We have marked the two points where this happens in Fig. [8] The mean 
bulk velocity is measured at a slightly different angle i/' + = 81 o ±9° and tf/ + (p = 261° + 8° 
with a bulk velocity of (280+7) km s _1 . This means t hat the offset angle is = -9° ± 7°. This 
position of the major axis agrees well with that found bv lSugerman et alj ( 12005b (PA= 8 1 .° 1 ± 0.°8) 
using imaging of light echoes and modelling the shape of the all the matter around the supernova. 
Our determination is not weigthed by intensity and rather than having single points we make a 
dynamical map of the whole ring. The influence of the different hot spots or the contamination 
by Star 4, which is directly superposed onto the inner stellar ring, appear negligible. We have 
marked the location of Spot 1 and Star 4 in the diagram and there are no clear deviations from the 
overall behaviour at these two poi nts. The derived velocity of the ring of (280+7) km s" 1 agrees 
with other results of 286.5 km/s ( Groningsson et alj|2007 
Cumming & Meiklelll993 ). 



Meaburn et al 



1995 



Crotts et al. 



1995 



Another result concerns velocity of the shocked ring material after passage of the shock. As we 
discussed above, the peak velocity is a function of the geometry of each shocked spot, as well as 
the direction of the ring plane to the line of sight. We can, however, use the geometric information 
above to find the deprojected, average shock velocity around the ring. Assuming an inclination 
angle of i = -42°. 8, and a systemic velocity of 286.5 km s _1 we have calculated the expansion 
velocity and binned the datapoints in Fig. [9] We derive from this a mean velocity around 90 km s _1 . 

Figure [10] shows the flux of the [Fen] A 1.257 //m, the Hei, and the Bry lines plotted vs. the 
azimuthal angle. We see here that the H and He lines, as expected, follow each other, and that 
the ratio He/H is fairly uniform around the ring. While the narrow (intermediate) [Fe n] line flux 
follows the He i flux level between 100° and 200° its flux increase around Spot 1 peaks at a different 
angle. This difference can be understood from the fact that while the H and He lines arise in an 
ionised region at ~ 10,000 K, the [Fe II] lines arise at ~ 5000 K, behind the ionized zone. The 
extent of this photoionized zone, and therefore the [Fe II] flux, depends on the shock velocity. In 
addition, collisional de-excitation is important for the [Fe II] lines, while the H and He lines are 
less sensitive to the density, being dominated by recombination. The ratio of the [Fe II] and H and 
He lines will therefore depend on both the shock velocity and density, which are both likely to vary 
along the ring. Because V s ~ V*biast(Po/Pbiob) the shock velocity into the cloud depends on the 
cloud density, pbiob- Finally, if the ratio of the narrow and intermediate components vary between 
the different lines this may result in different peak velocities. 

We expect the flux of H and He to increase around the ring as the hot spots evolve. 

We will determine the evolution of the fluxes, as well as the velocity, around the ring in our 
next paper with epoch 2005 SINFONI data. 



4. Summary 

The supernova shock interaction with the inner circumstellar ring around SN 1987A is in full swing. 
Our integral-field observations have revealed the extent of individual emission lines around the ring. 
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Fig. 9. Deprojected velocity vs. angle for PayS, Bry and He i. degrees is North, the angle increases 
from North to East. The curves represents the propogated uncertainties for the deprojection, and 
we have cut the data where the uncertainty was larger than the velocity measurement. In the depro- 
jection we assumed a systemic velocity of 286.5 km/s. 



From the high-resolution optical spectr oscopy of the coronal lin es from UVES dGroningsson et al 



2006) and the soft X-ray component (IZhekov et a l. 2005, 2006), we know that the shock is heating 



the ring material to very high temperatures. The material behind the shock is cooling rapidly, and 
the temporal evolution of the lines will tell us more about the conditions in these post shock areas. 
Our SINFONI Science Verification data only show the strongest lines and we are now investigating 
deeper integrations from a year later. The comparison of individual lines will be reserved for a later 
paper. 

Our data are still consistent with a circular ring, and we have derived an average recession 
velocity, which agrees with other measurements. Yet, we find that there are significant differences 
of the velocities around the ring, with the eastern part showing larger velocities than the western 
segment. T his corresponds well w ith the increase in surfa ce brightness see n in both the radio ob- 



Spot 1 in the East dMichael et al 



servations (Manchester et al. 2002 ), Chandra observation s (Park et al. 2006) and the appearance of 



2002; 



Pun et al 



2002; 



Sugermanetal 



20021) . By now the ring 



has been lit up in spots rather evenly spaced around the ring, which we also observe in the emission 
lines of H and He. 

From the velocity shift of the lines we find an average velocity of ~ 100 km s _1 for the bulk 
of the H and He emitting shocked gas. We emphasise that this is just the average velocity, and 
considerably higher velocity sh ocks are present, as is directly evidenced by the UVES observations 



dGroningsson e t al. 2006, 



2007). 




Fig. 10. Flux variations for different angular positions around the ring. The aperture area is A=n * 
(0.375") 2 for the H and He lines, and A=tt * (0.35") 2 for the [Fe n] line. 



The emission from the supernova itself has faded continuously. Today it is mostly the radioac- 
tive decay of 44 Ti and reverse shock interaction, which powers the emission from the ejecta. With 
the integral-field spectrograph we can actually spatially isolate the ejecta emission from the ring 
emission. The only clear line we detect from the ejecta is [Fen] 1.644 jjm with a width of about 
3400 km s _1 . This is consistent with the radioactively and reve rse shock heated ejecta seen in 



Ha, Ca II and Mg I (S mith et al 



2005 



Groningsson et al.ll2007l) . Deeper observations are likely 



to detect more lines, as well as the broad wings from the reverse shock in the H and He lines. An 
advantage of the IR lines is that the [Fe n] lines are less blended than in the optical spectrum. 

Monitoring the transformation of SN 1987 A into a supernova remnant is a unique opportunity. 
It is the only such object in which we can directly observe this event while simultaneously resolving 
the interaction region. The structure of the circumstellar ring has been a puzzle since its discovery 
when it was ionised by the SN soft X-ray flash. For a complete picture we need to monitor this 
transition at all wavelengths. The SINFONI observations we are obtaining are unique in that they 
allow us to spatially investigate the emission sites around the ring and infer the local conditions of 
the interaction. 
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